Residual absorption at zero temperature in d-wave superconductors 
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In a d-wave superconductor with elastic impurity scattering, not all the available optical spectral 
weight goes into the condensate at zero temperature, and this leads to residual absorption. We find 
that for a range of impurity parameters in the intermediate coupling regime between Born (weak) 
and unitary (strong) limit, significant oscillator strength remains which exhibits a cusp like behavior 
of the real part of the optical conductivity with upward curvature as a function of frequency, as 
well as a quasilinear temperature dependence of the superfiuid density. The calculations offer an 
explanation of recent data on ortho-II YBa2Cu306.s which has been considered anomalous. 
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I. INTRODUCTION 



Recent measurements^ of the microwave conductivity up to 21 GHz at low temperatures in single crystal oxygen 
ordered ortho-II YBa2Cu30g.5 (YBCO) have revealed that there remains significant absorption at T — in the 
superconducting state, while at the same time the inverse square of the penetration depth shows a linear in temperature 
behavior. In the regime where this happens, the real part of the optical conductivity ai(u>) shows an unusual 
characteristic upward curvature very different from a Drudedike form. Anomalous residual absorption has also been 
seen in the past in THz work on thin films of Bi2Sr2CaCu20s+5 (Bi2212)2 and in this case was interpreted as due 
to coupling to a collective mode, with the spectral weight of this mode tracking the temperature dependence of the 
superfiuid density. 

In a pure <i-wave BCS superconductor all of the optical spectral weight will go into the condensate at T = providing 
a delta function to the real part of the conductivity at zero frequency. This follows from the Ferrell-Glover-Tinkham 
sum rule2i4 which holds when the kinetic energy change between normal and superconducting state is negligible. 
When elastic impurity scattering is introduced into the system, this no longer holds. Some optical spectral weight 
will remain at finite frequency under the real part of the conductivity at T — and correspondingly, the superfiuid 
density will be reduced below its pure limit value. At the same time, however, in the unitary (strong) limit of 
impurity scattering which is the most studied case so far, a switch from linear to quadratic behavior in temperature 
is expected for the penetration depth. In this work we consider the case when the impurity scattering strength falls 
in the intermediate regime between Born (weak) and unitary (strong) limit. We find that for a range of impurity 
parameters, i.e.: concentration and strength of the potential measured in terms of the width of the electronic band, 
the residual, impurity induced absorption remains very significant, while at the same time the superfiuid density obeys 
a quasilinear law over the entire low temperature range and, more importantly, crj (to) shows anomalous behavior in 
agreement with the experimental observation in ortho-II YBCO. 

II. FORMALISM 

The general expression for the optical conductivity at temperature T and frequency v is well known and has been 
given in many places It can be written as: 



oc oo 

— J did tanh Je{<^, v) + J did tanh 

-v 



0(u + v) 



Je(-cJ - v, v) 



(1) 



where Q p is the plasma frequency and the function Jg(cd, v) takes on the form 



2J e (w,v) = 1 [l-N(9,u)N(0,u + v)-P(9,u)P(9,u + u)] 
hj\ + hj-2 

1 [1 + N*(9,ui)N(9, oj + v)+ P*(9,ui)P(6,ui + v)] . (2) 
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In Eq. (3 — 1/ksT, with Ub the Boltzmann factor and T the temperature. In Eq. (J2J 

E x {9,uj) = ^Co 2 (id + iO+)- A 2 (6»,w), 
E 2 {9,w,v) = E 1 (9,u> + v), 

and 

N(9^) = ^±^l, P(9,.)- K ^ 



Here 9 is the angle giving the direction of momentum on the Fermi surface in the two dimensional Cu02 Brillouin zone. 
In BCS the gap does not depend on frequency u> and has the form A(9,oj) = A(T) cos(2#) with A(T) a temperature 
dependent amplitude. The angular bracket (• • • )g indicates an average over angle 9. 

The renormalized frequencies ui(a}) do not depend on angle but are changed because of the elastic impurity scattering 
which is the only scattering processes accounted for. The equation for o)(w+iO + ) in a i-matrix approximation igSiSiSiifl 

u>(u> + z0+) = u, + »^r+ {N ^\ 2 . (3) 
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r + is proportional to the impurity concentration and c = 1/ (2irN(0)Vi mp ) with N(0) the normal band density of 
states at the Fermi surface and Vi mp the impurity potential. An important parameter is 7 defined as 3m(i(w + i0 + ) 
at us = which is shown in the bottom frame of Fig. 01 For small values of Vi mp , Born's approximation applies and 
c — > 00 with t + = T + /c 2 and the impurity term in Eq. becomes proportional to (N(9,us))g, its real part is the 
superconducting state density of states (N s (us)), while for large values of Vi mp , (c — > 0) it is inversely proportional 
to (N(9, us)) s ) (unitary limit). Both Born and unitary limit are extreme and a finite c is more realistic. As N(0) 
is roughly inversely proportional to the electronic band width (W), we have Vi mp ~ W/(2irc). For c = 0.2 the 
corresponding Vi mp is 1.3 times the band width. 



III. RESULTS AND DISCUSSION 



The quasiparticle density of states N s (u>) for a <i-wave superconductor in the case of a finite value of c — 0.2 is 
shown in Fig. ^ (top frame). The gap amplitude was set at 24v / 2meV and the impurity parameter r + was varied 
from T+ = 0.01 mcV (solid), T+ = 0.05 meV (dashed), T+ = 0.1 meV (dotted), and T+ = 0.15 meV (dash-dotted). 
Clearly, impurities affect most strongly the low energy region near us — up to, say, 10 meV. The pure limit, a straight 
line through the origin is approached as T + is reduced, but even for the smallest value of T + = 0.01 meV, changes in 
N s (us) persist up to us ~ 4 meV. It is the behavior of N s (us) in the low energy range which determines low temperature 
properties such as the specific heat. The behavior found for a finite c could not be inferred from a knowledge of the 
unitary or Born limit. In the Born limit the deviation from linearity would only occur at an exponentially small value 
of us. Results for the unitary limit are shown in the middle frame and are seen to be very different from the top frame. 
The energy dependence of the quasiparticle lifetime of Eq. © 3mw(w + i0 + ) is shown in the bottom frame of Fig. ^ 
For c = 0, the scattering rate has its maximum at us — and then drops monotonically as u> increases. For the larger 
values of c shown here this is no longer the case. A maximum in scattering rate occurs around us ~ 8 meV which is 
much larger than its value at us = (denoted 7). 

The real part of the optical conductivity for a <i-wave BCS superconductor is shown in Fig. [3 in units of Slp/87r. 
Five values of temperature are considered, namely T = OK (short dotted line) T — 5K (solid), T = 10 K (dashed), 
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2: The real part of the optical conductivity <ti(T, lo) vs uj for five temperatures. The gap amplitude Ao = 24\/2meV and 
the impurity parameter F + = 0.01 meV. The top frame is for c = 0.3 and the bottom for c = 0. The computer units on the 
conductivity are such that a factor of f2p/87r needs to be applied. 



FIG. 



T = 15 K (dotted), and T — 20 K (dash-dotted). The d-wave gap amplitude was set at 24v / 2meV. The top frame 
applies to an intermediate scattering case with T + = 0.01 meV and c = 0.3. The bottom frame is for comparison and 
shows unitary limit results (c = 0). The real part of the optical conductivity a\{T,uj) vs u> up to lmeV, is radically 
different in the two cases. For c = 0.3 the finite T curves are concave up and drop rapidly with increasing ui and show 
a very narrow half-width in qualitative agreement with the experiments by Turner et alX and Hosseini et alil> The 
concave upward nature of these curves is due to the energy dependence of the quasiparticle scattering rate obtained 
from Eq. J3J). These curves all correspond to the temperature dominated regime T>7 (see Fig. [5J bottom frame). 
Their half width is not related to the value of 7 but rather to temperature which samples the important energy 
dependence in Smiifw + i0 + ). The T = curve is impurity dominated and, consequently, shows a qualitatively 
different behavior. It cannot be obtained from an extrapolation to T = of the finite T curves shown. These 
conductivity curves are in qualitative agreement with the experimental findings^ In the unitary limit the curves show 
instead a concave downward behavior at small ui and a broad half-width. This can be traced to the fact that the 
zero frequency scattering rate 7, is a rapidly increasing function of decreasing c as is seen in the bottom frame of 
Fig. [5| For the unitary limit 7 is greater than its normal state value 7rr + and is given by the well known relation 
7 ~ 0.63vwAoT+ while in the Born limit 7 is exponentially dependent on T + (7(0 — > 00) = 4Ao exp(— Aq/2Tn) with 

r w = r+/c 2 .) 

The quantity of interest is the total spectral weight remaining under the real part of the optical conductivity at 
temperature T, A(T) — f^div <7i(T, ui). Numerical results are presented in Fig. [3] which consists of three frames. The 
top frame is for c = 0.3 while the middle frame, included for comparison, is for the unitary limit. In both frames 
A = 24v / 2meV and the temperature range considered is up to 20 K (t = T/T c ~ 0.11). The solid squares are 
for r + = 0.15 meV, solid circles T + = 0.1 meV, and solid triangles T + = 0.05 meV. The curves show a quasi linear 
temperature dependence over the entire temperature range considered. This is in sharp contrast to the results in the 
middle frame. The temperature dependence is now quadratic which is the behavior expected in the unitary limit. In 
the bottom frame we compare A(T) vs T with the corresponding results for the penetration depth A~ 2 (0) — A _2 (T) 
vs T. We see that the two curves are parallel and that the dotted curve goes through zero at the origin. This 
follows directly from the Ferrell- Glover- Tinkham sum rule which is satisfied to within our numerical accuracy. In our 
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FIG. 3: Top frame: the remaining spectral weight A(T) — f Q °^ du <ti (T, uj) as a function of temperature T in the range to 20 K 
for gap Ao = 24\/2meV. The impurity potential is characterized by c — 0.3 (to be compared with c = for the unitary limit 
in the middle frame). The bottom frame compares (2/w)A(T) with A" 2 (0) - A _2 (T) for the case F =0.1 meV and c = 0.3. 



computer units the sum rule is 

oo 

lim ua 2 (T, u) + - [ du a x (T, w) = 2, (4) 

0+ 

where 02 is the imaginary part of the conductivity. While we have not attempted a best fit to experiment, the results 
of the bottom frame of Fig. [3] are in qualitative agreement with the findings of Turner et a/i 

Zero temperature values for A(Q) are shown in the top frame of Fig. 0] The lines are a guide to the eye through 
numerical results based on full conductivity calculations. The points are solid down-triangles for T + = 0.3 meV, 
solid squares for T + — 0.15 meV, solid circles for T + = 0.1 meV, solid up-triangles for T + = 0.05 meV, and solid 
diamonds for T + = 0.01 meV. In all cases A(0), which is entirely due to impurities, drops with increasing c. It is 
largest in the unitary limit but remains significant even when c has increased to 0.4 provided T + is not too small. In 
the bottom frame we show 7(c) vs c, and stress that it becomes small more rapidly than does A(0) as c increases; 
consequently, significant absorption can remain even with small 7. As can be seen in the top frame of Fig. [2 the 
value of A(0), which gives the optical weight under <j\(u>), samples importantly the quasiparticle scattering rate away 
from lu — and, thus, A(0) is not directly related to 7(c). There is much more absorption than would be indicated 
by the value of 7(c) at the larger values of c considered here. This is emphasized by the solid gray line in the top 
frame of Fig. 21 which corresponds to the results of the approximate formula ^4(0) ~ (7/A0) ln(4Ao/7) obtained by 
approximating the w-dependent quasiparticle scattering rate in the exact expression by its lu — value, i.e.: by 7(c). 
Here T + = 0.15 meV. Comparison with the full results (solid squares) reveals that the agreement is reasonable for 
c < 0.1, i.e.: the region for which the quasiparticle scattering rate has its maximum at 10 — 0, but fails when its 
maximum falls instead at finite 10 (see Fig.^ bottom frame). 

In their analysis of optical conductivity of thin films of Bi2212 Corson et al. 2 found a remaining spectral weight at 
T = of about 30% of the corresponding value of the superfluid density. In our language this corresponds to a value 
of ,4(0) ~ 0.5. Corson et al. 2 also found that the remaining optical weight was distributed over a width of order 12 K 
independent of temperature. Reference to the theoretical curves in Fig. 0] shows that this value of ,4(0) (top frame) 
is rather high for the observed value of the width ~ 7 (bottom frame); the data suggests c ~ 0.4 and T + of the same 
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FIG. 4: Top frame: remaining spectral weight A(0) = duo g\(Q,lu) as a function of the impurity potential strength c. The 
bottom frame gives the zero frequency value of the effective scattering in the superconducting state 7(c) as a function of c. The 
solid gray line represents the results of the approximate formula for Y + — 0.15 meV. 

order as was required^ to fit the single crystal infrared optical data of Tu et alH in Bi2212, T + = 0.61 meV. While 
Corson et al. interpret their data in terms of extra absorption due to a collective mode, part if not all is due instead 
to the impurity induced absorption discussed in this paper. 

We mention previous work on film of YBa2Cu30y_5 by Hensen et al. 1A in which the sizable surface resistance 
observed at low temperatures and attendant quasi linear temperature dependent superfluid density is seen to be 
consistent with <i-wave and impurity scattering with a phase shift intermediate between Born and unitary limit. This 
work, however, is limited to the consideration of a single microwave frequency, namely 87 GHz, and thus remains 
silent on the unusual frequency dependence of a\ (lu) vs lo found also to be the characteristic of this regime in our 
theoretical work and in the experimental data of Turner et alX 

IV. CONCLUSION 

For impurity scattering with a phase shift intermediate between and 7r/2 we have found a regime in which consid- 
erable residual absorption can remain at zero temperature while at the same time the superfluid density is quasilincar 
in temperature. In addition, and more importantly, in this regime, the real part of the optical conductivity, <j\(lo), 
exhibits a distinctly non-Drude dependence on frequency in the microwave region. At low, but finite, temperatures 
<J\{uj) vs to has a cusp like dependence on lo out of lu = 0, followed by a rapid drop (concave upward curvature) 
as lo increases. The width at half maximum of ai(u)), which can be very small, varies strongly with temperature. 
This behavior is traced to the strong energy dependence of quasiparticle lifetime which is sampled in the temperature 
dominated regime discussed in this work. All aspects of this work are in good qualitative agreement with recent data 
in ortho-II YBa2Cu306.5- 
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